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We present a dynamical description and analysis of non-equilibrium transitions in the noisy one- 
dimensional Ginzburg-Landau equation for an extensive system based on a weak noise canonical 
phase space formulation of the Freidlin-Wentzel or Martin-Siggia-Rose methods. We derive prop- 
agating nonlinear domain wall or soliton solutions of the resulting canonical field equations with 
superimposed diffusive modes. The transition pathways are characterized by the nucleations and 
subsequent propagation of domain walls. We discuss the general switching scenario in terms of a 
dilute gas of propagating domain walls and evaluate the Arrhenius factor in terms of the associated 
action. We find excellent agreement with recent numerical optimization studies. 
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I. INTRODUCTION 



Phenomena far from equilibrium are ubiquitous, including turbulence in fluids, interface and growth problems, 
chemical and biological systems, and problems in material science and nanophysics. In this context the dynamics of 
' complex systems driven by weak noise, corresponding to rare events, is of particular interest. The issue of different 
t^J- , time scales does in fact characterize many interesting processes in nature. For instance, in the case of chemical 
reactions the reaction time is often orders of magnitude larger than the molecular vibration periods Q. The time 
scale separation problem is also encountered in the case of conformational changes of biomolecules, nucieation events 
during phase transitions, switching of the magnetization in magnetic materials 0! Q , and even in the case of comets 
exhibiting rapid transitions between heliocentric orbits around Jupiter Q. 

The weak noise limit, characterizing the time scale separation, is associated with the strong coupling regime, and 
the problem of determining kinetic pathways and transition probabilities between metastable states in systems with 
many degrees of freedom presents one of the most important and challenging tasks in many areas of physics [5| . The 
long time scales are associated with the separation in energy scales of the thermal energy and the energy barriers 
between metastable states; the transition takes place by sudden jumps between metastable states followed by long 
waiting times in the vicinity of the states. The fundamental issue is thus the determination of transition pathways 
and the associated transition rates. 

In the weak noise limit the standard Monte Carlo method or direct simulation of the Langevin equation becomes 
impractical owing to the large separation of time scales and alternative methods have been developed. The most 
notable analytical approach is the formulation due to Freidlin and Wentzel which yields the transition probabilities 
in terms of an action functional This approach is the analogue of the variational principle proposed by Machlup 
and Onsager 0, Q, see also work by Graham et al. 0, [ljj. The Freidlin-Wentzel approach is equivalent to the 



iproacn is equivaiei 

Martin-Siggia-Rose method [TJJ in the weak noise limit of the path integral formulation Il3l HH \wi fl6| . In order 
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to overcome the time scale gap various numerical methods have also been proposed. We mention here the transition 
path sampling method 01 an d optimization techniques [l8l . 

A particularly interesting non-equilibrium problem of relevance in the nanophvsics of switches is the influence of 
thermal noise on two-level systems with spatial degrees of freedom, see Refs. [2 EES • In a recent paper by E, Ren, 
and Vanden-Eijden [2l|, see also Ref. this problem has been addressed using the Ginzburg-Landau equation 

driven by thermal noise. Applying the field theoretic version of the Onsager-Machlup functional m the Freidlin- 
Wentzell formulation Q , these authors implement a powerful numerical optimization techniques for the determination 
of the space-time configuration minimizing the Freidlin-Wentzell action and in this way determine the kinetic pathways 
and their associated action, yielding the switching probabilities in the long time-low temperature limit. The picture 
that emerges from this numerical study is that of noise-induced nucleation and subsequent propagation of domain 
walls across the sample, giving rise to the switch between metastable states. 

In recent work we have addressed a related problem in non-equilibrium physics, namely the Kardar-Parisi-Zhang 
equation or equivalent noisy Burgers equation describing, for example, a growing interface in a random environment. 
Using a canonical phase space method derived from the weak noise limit of the Martin-Siggia-Rose functional or 
directly from the Fokker-Planck equation, we have, in the one dimensional case, analyzed the coupled field equations 
minimizing the action both analytically |22l l23l [2J, l25l 1261 127| and numerically [23] • The picture that emerges is that 
the transition probabilities in the weak noise limit are associated with soliton propagation and nucleation resulting 
from soliton collisions. 

In the present paper we apply the canonical phase space approach to the noisy Ginzburg-Landau equation discussed 
by E, Ren, and Vanden-Eijden |2l| an d attempt to account for some of their numerical findings. We thus give 
analytical arguments for the propagation of noise-induced solitons and the nucleation events originating from soliton 
annihilation and creation. The paper is organized in the following way. In Sec. II we introduce the noisy Ginzburg- 
Landau equation. In Sec. Ill we review the canonical phase space approach. In Sec. IV we discuss diffusive mode 
and domain wall solutions of the field equations replacing the noisy Ginzburg-Landau equation. In Sec. V we analyze 
the domain wall dynamics. In Sec. VI we present a stochastic interpretation of our results. In Sec. VII we discuss 
equilibrium properties and kinetic transitions. In Sec. VIII we present an interpretation of the numerical results 
obtained by the optimization studies of E et. al 1211 . Sec. IX is devoted to a summary and a conclusion. A brief 
version of the present work has appeared in Ref. [2{j . 



II. THE NOISY GINZBURG - LANDAU MODEL 



The noisy Ginzburg-Landau equation driven by white noise has the form 

^ = -r— + J?, 2.1 

at ou 

where the locally correlated Gaussian white noise is specified by the moment 

( V (x,t)r,(0,0)) = AS(x)S(t), (2.2) 
with noise strength parameter A. The free energy F providing the deterministic drive is given by 

'-5/*((*)'-H- (23) 

In the switching problem considered by E et al. [2lJ the double well potential V(u) has the form 

V(u)=k%(l-u(xf) 2 , (2.4) 

with strength parameter k^. The potential vanishes at the two minima u ± 1 and assumes the maximum value 
V(0) = /cq at u = 0. The time scale is set by the kinetic transport coefficient F and the dimensionless scalar field 
u determines the spatial and temporal state of the switch. The "Mexican hat" double well potential is depicted in 
Fig.ffl 

The Ginzburg-Landau equation in its deterministic form has been used both in the context of phase ordering 
kinetics [30j and in its complex form in the study of pattern formation |3l| . In the noisy case for a finite system the 
equation has been studied in j^]; see also an analysis of the related </> 4 theory in [33J. In the present problem the 
noisy equation provides a generalization of the classical Kramers problem 5] to spatially extended systems. 

By inspection of Eqs. I]2.ip . (|2.2|) . and (|2.3|l we note that ko has the dimension of an inverse length, T the dimension 
of length squared over time, and A the dimension of velocity. For large ko, which is the case considered here, and 
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FIG. 1: We depict the double-well potential defining the unperturbed states of the switch. The potential has minima at u — ±1 
and a maximum at u — 0. 



for ?y = the potential term dominates the diffusive term and the field u locks on to the values ±1 in bulk; imposed 
boundary values being accommodated over a saturation or healing length of order A;^" 1 . In the presence of noise the 
minima become unstable to thermal fluctuations and the noise can drive the system over the potential barrier. 
The Ginzburg-Landau equation 1|2.1[1 admits a fluctuation-dissipation theorem yielding the stationary distribution 



P sta t oc exp 



(2.5) 



In a thermal environment at temperature T we have A = 2TT (in units such that fee = 1) and P s tat cx exp[—F/T], 
i.e., the Boltzmann distribution. The equilibrium states follow from the condition 

SF 

— = -V 2 u- 2fc 2 u(l - u 2 ) = 0, (2.6) 
ou 

yielding two degenerate uniform ground states u = ±1 with free energy F = 0, as well as two nonuniform domain 
wall solutions 

u dw (x) = ±taiihk (x - x ), (2.7) 

centered at Xq and of width fc^" 1 . For large \x\ the domain wall solutions overlap with the uniform ground state 
solutions u — ±1. The associated domain wall free energy is 

F dw = 4fc /3. (2.8) 

The minima, maxima, and saddle point structure of the complex energy landscape of F as a functional of {u(x)} are 
inferred from the spectrum of the differential operator 

^1 = -V 2 -2fc 2 (l-3u 2 ). (2.9) 

For the ground states u = ±1 the spectrum of 5 2 F/Su 2 is given by the plane wave mode ~ exp(ikx) with positive 
eigenvalues k 2 + 4fcg. For a single domain wall configuration or n connected domain walls the spectrum is com- 
posed of zero-eigenvalue translation modes (Goldstone modes) and superimposed phase shifted plane wave modes 
~ exp[ikx) exp(i(j)) with positive eigenvalues k 2 + 4fcg, see e.g. Ref. We thus infer that the free energy landscape 
possesses two global minima at u = ±1 and a series of local saddle points of free energy 4nfco/3, corresponding to n 
connected domain walls at positions Xi,i = 1, 2, ■ • ■ n. In Fig. [3 we have depicted a static 3-domain wall configuration 
connecting the ground states u = — 1 and u = +1 with free energy 4fco. The static single domain wall or multi-domain 
wall configurations correspond to points in the free energy landscape F = F({u(x)}) and cannot effectuate a switch 
between the ground state configurations u = ±1. However, in the presence of noise the Kramers escape mechanism 
sets in and the system can perform a kinetic transition. The transition probability is typically characterized by an 
Arrhenius factor P ~ exp[— TSF/A] where SF is the free energy associated with the potential barriers encountered in 
the dynamic transition, see e.g. In order to address this issue for the noisy Ginzburg-Landau equation we apply 
the canonical phase space approach. 
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FIG. 2: We depict a static 3-domain wall configuration connecting the ground states u = —1 and u — +1. The domain walls 
are located at x\, X2, and £3. The configuration corresponds to a saddle point in the free energy landscape and has free energy 
4fc . 



III. THE CANONICAL PHASE SPACE APPROACH 



The canonical phase space approach is discussed in detail in |28l l35j . Briefly, the Fokker-Planck equation for the 
probability distribution P({u(x)},t) has the general form 



dP 



(3.1) 



with formal solution P oc exp[fft/A]. The Fokker-Planck equation 1|3.1[1 is driven by the Hamiltonian or Liouvil- 
lian H(u,S/Su) (a differential-functional operator in u(x) space). The method is based on a weak noise WKB-like 
approximation 



P oc exp 



S_ 
A 



(3.2) 



applied to Eq. (|3.1|l . To leading order in the noise strength A the action S satisfies the Hamilton- Jacobi equation 

(3.3) 



— +H(p,u) = Q, 



with canonical momentum p — SS/du, defining the Hamiltonian H as a functional of u and p and implying a principle 
of least action. The variational principle 5S — 0, moreover, implies the Hamiltonian field equations of motion 



du _ SH 

'm^lp' 1 

dp _ 6H 
dt Su ' 



and the action 



S = I dxdt p 



du 
dt 



dt H. 



Applying this scheme to the Ginzburg-Landau equation 12. 1|) we obtain 



du 

~dt 



6F 

-r— + P , 

ou 



dp S 2 F 
dt ou z 



driven by the Hamiltonian (the generator of time translations) 



H = J dxH=^J dxp 



p-2T 



Su 



(3.4) 
(3.5) 

(3.6) 

(3.7) 
(3.8) 

(3.9) 



5 



with SF/du and S 2 F/8u 2 given by Eqs. 1)2.6)1 and (|2.9|l . As a result the action S associated with an orbit from u\ to 
«2 traversed in time T is given by 



ru 2 ,T 

S(ui,U2,T)= / dxdt 



du 



(3.10) 



where TL is the Hamiltonian density. In our further analysis we shall also make use of the total momentum II (the 
generator of space translations) 

11= f dx u^r-. (3.11) 
J dx 

The prescription for determining the relevant Arrhenius factor associated with a transition and the kinetic pathway 
is thus straightforward. Fixing boundary values in space and time the first step is to solve the coupled field equations 
(|3.7() and 1)3.8)1 for an orbit (a minimizer in the terminology of Ref. from u\ to U2 traversed in time T, next we 

evaluate the action S along an orbit according to Eq. I|3.10|l . and finally deduce the transition probability from Eq. 

The correspondence with the equivalent Freidlin-Wentzel action Sfw is obtained by inserting the equation of motion 
(|3~7|) and the Hamiltonian (|3~3|) in the action (|3~TU|) . yielding S = (1/2) J dxdtp 2 or S F w = (1/2) Jdxdt(du/dt + 
TSF/Su) 2 . The Freidlin-Wentzel method is a Lagrangian configuration space method with Lagrangian density C = 
(l/2)(du/dt + TSF/6u) 2 , whereas the present approach is a phase space formulation with action l|3.10|l . We also note 
that the equations of motion 1)3.7)) and 1)3. 8|l are identical to the saddle point equations in the Martin- Siggia-Rose 
functional formulation 0, . 

IV. DIFFUSIVE MODE AND DOMAIN WALL SOLUTIONS 

More explicitly, inserting SF/Su and 5 2 F/Su 2 from Eqs. 1)2.6)1 and (|2.9|l the equations of motion l|3.7|) and l|3.8|) 

assume the form 

^ = rg + 2Ffc Vi-- 2 )+P, (4.1) 
(>p = -r§-2rfc 2 Ki-3u 2 ), (4.2) 



dt dx 2 



and the Hamiltonian is 



H=\J ' dxp(p + 2r0 + 4Tk 2 u(l - u 2 )^j . 



(4.3) 



The equations of motion (|4.1(l and (|4.2(l determine orbits in a multi-dimensional phase space {u(x),p(x)} lying on 
the energy manifolds determined by Eq. (|4.3() : for open or periodic boundary conditions the orbits are moreover 
confined by the conservation of the total momentum n given by Eq. (|3.11(l . From an analytical point of view the 
field equations 1|4.1JI and 1)4. 2JI are in general intractable. Even numerically, the negative diffusion term in Eq. 14.2)) 
renders the coupled equations highly unstable as was noted in the numerical analysis of the Burgers equation |2a | . 

A. Linear diffusive modes 

It is instructive first to consider the easily discussed linear case of simple diffusion for small k$. Thus ignoring the 
non-linear potential terms in Eqs. ()4.1)) and ()4.2)) we obtain the linear equations 

du ^d 2 u 

m- T dx- 2+ ^ (44) 

dp d 2 p 



generated by the Hamiltonian 



H = ~ / il.rp 



dx 2 



(4.6) 
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P k =2rk 2 u k 




FIG. 3: We show the orbits in the (uk,Pk) phase space. The finite time orbit from un, to «2fc lies on the energy manifold 
Ek = ptlPk — 2Yk 2 Uk\- In the long time limit T — > oo the orbits migrate to the transient manifold pk = (I) and the stationary 
manifold pk = 2Vk 2 Uk (II) passing through the saddle point (iik,Pk) = (0,0) (SP) implying ergodicity and a stationary state. 



In Fourier space, setting Uk — J dx exp (-ikx)u and pk = J dx exp(—ikx)p, u* k = u- k and p k = p-k, the equations 
of motion decompose and we arrive at the orbit solution, see Ref. |3jj | ■ 

u 2k sinh Tk 2 t + u lk sinh Tk 2 (T -t) 



u k (t) 



sinhrfc 2 T 



(4.7) 



from uik at time t = Q to u 2 k at time t = T; u% t 2 — j{dk/2ir) e-xjp(ikx)ui.2, k - The noise field pk is slaved to the motion 
of Uk and given by 



-Tk A T 



Pk (t) = rfcV 



sinhrfc 2 T 

Likewise, the Hamiltonian H decomposes into independent fc-mode contributions 

H = j ^p* k [pk-2Tk 2 u k }. 

Inserting the solutions l|4.7|l and l|4.8|l we obtain specifically for the energy of the fc-th mode 

(rfc 2 ) 2 \u 2k \ 2 + \u lk \ 2 - 2u lk u 2k coshrfc 2 T 



Ek 



sinh 2 rfc 2 T 



(4.8) 



(4.9) 



(4.10) 



The orbits lie on the energy manifolds given by Ek = PkiPk — 2Tk 2 Uk\- In the long time limit T — > oo the orbits 
migrate to the zero-energy manifolds consisting of the transient submanifold pk = and the stationary submanifold 
Pk = 2Tk 2 Uk, and pass asymptotically through the hyperbolic fixed point at (uk,Pk) — (0,0), determining the 
stationary state. In Fig. [3] we have shown the orbits for a particular fc-mode in a plot of p k versus u k . 
Finally, the conserved momentum II and the action 5* follow from Eqs. I|3.11fl and i|3.10|) 



n = J ^i~ ik )p*k u k, 



and 



<* rfc2 l u 2fc - u lk e 
2ir 



-Tk £ T\2 



1 



2Tk 2 T 



yielding the transition probability from u\k to u 2k in time T 

P(n lh . no,, . T) x exp - — 



dk 2 \u 2 k - uike' 



-Tk 2 T\2 



(4.11) 



(4.12) 



(4.13) 
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In the long time limit T 



oo we reach the stationary distribution 

1 f dk 



P s t({u k }) oc exp 



2tt 



Tk A \u k \ 



(4.14) 



This result is consistent with the noise driven diffusion equation written in the form du/dt = —TSF/Su + r), (7777) = 
AS(x)5(t), with free energy F — (1/2) J dx(du / dx) 2 . Invoking the fluctuation-dissipation theorem we then obtain an 
equilibrium distribution given by the Boltzmann factor exp(— F/T), A = 2TT, consistent with Eq. (|4.14ll . 

In summary, in the linear purely diffusive case the configurations u(x, t) decompose into independent fc-modcs. 
In the dynamical phase space approach the noise r\ k driving the individual /c-modes is replaced by the noise field 
Pk which couples parametrically to the time evolution of the u k mode. In a transient time regime u k is damped 
according to the diffusion equation and the orbit in phase space lies close to the transient or noiseless submanifold 
Pk = and approaches the saddle point (uk,Pk) = (0,0). At longer times the growing noise field pk drives u k away 
from the saddle point and the orbit approaches the stationary or noisy submanifold p k = 2Tk 2 Uk, i.e., the distribution 
associated with the noise-driven diffusion equation approaches a stationary distribution. In the limit T — > 00 the 
orbit passes asymptotically through the saddle point and the orbit from tin to 11,2k lies on the zero-energy manifolds 
determining the stationary state. 



B. Domain wall modes 



In the nonlinear case the phase space representation of the noise-driven Ginzburg Landau equation is given by the 
coupled field equations i|4.1|) and Ij4.2|l determining the orbits. As in the linear case we can identify the zero-energy 
submanifolds determining the stationary state. This is related to the existence of a fluctuation-dissipation theorem 
for the noisy Ginzburg-Landau equation (|2.1|l expressed by the existence of a free energy. From the Hamiltonian 
(|3.9() we infer tentatively the zero-energy submanifolds p = and p = 2T5F/Su. The transient noiseless submanifold 
p = is consistent with the equations of motion and the configurations decay according to the damped deterministic 
Ginzburg-Landau equation 

Note that in contrast to the linear case where the u configurations decompose in fc-modes which decay according to 
itfc oc exp(— Yk 2 t), an initial configuration u\ in the nonlinear case will in general decay forming a pattern of interacting 
and annihilating domain wall configurations with superimposed diffusive modes. 

The stationary submanifold p = 2TSF/8u inserted in Eqs. I|3.7[l and (|3.8|l yields the equations 

du SF , 

dp ^ 2 5 (SF 



«=^UJ' (4 - 17) 

which are consistent since du/dt — (l/2)p. It is interesting to notice that the motion of u on the stationary noisy 
submanifold p = 2T(5F/ Su) is a time-reversed version of the motion on the transient submanifold p = 0. Finally, on the 
zero-energy manifold the action in Eq. (|3.t)[) takes the form S = J dxdt pdu/dt = 2T J dxdt (du/ dt)(8F / Su) = 2TF, 
yielding the stationary distribution in Eq. I|2.5I) for the noisy Ginzburg-Landau equation. 

In the absence of a fluctuation-dissipation theorem as is for example the case for the kinetic KPZ equation or, 
equivalently, noisy Burgers equation, we can in general not explicitly identify the stationary zero-energy submanifold 
and thus simply determine the stationary state; an exception is the KPZ equation in ID where for special reasons a 
fluctuation-dissipation theorem is available, see Refs. |23l Il2 lL l35|. 

In order to address non-equilibrium properties such as a specific transition probability from an initial state ui to a 
final state 1x2 in passage time T we must address the nonlinear equations of motion jOJ and (|4.2|) . The field equations 
are not integrable and do not yield to a general analytical solution. However, we can advance our understanding by 
first searching for static solutions on the transient manifold p = 0, i.e., the solution of the equation SF/Su = yielding 
according to Eq. 12.6fl the domain wall solutions in Eq. 12.7fl . The domain wall excitations are of the instanton type 
and can be located at arbitrary positions, see Ref. j3g. Since the overlap between two well-separated domain walls 
is exponentially small we can construct approximate multi-domain wall solutions of the form 



= t7j tanhfco(x — Xj) + up, (4-18) 



i=l 



Pdw = 0. (4.19) 
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Here the parity index ui = ±1 for right hand and left hand domain walls, respectively, and x\ indicates the center of the 
domain wall. The offset u = for multi-domain wall configurations overlapping for large \x\ with two different ground 
state configurations and uq = ±1 for configurations overlapping with identical ground states u = ±1, respectively. 
Assuming that the inter-domain wall distance \xi — is large compared with the domain wall width 1/fco, i.e., the 
case of a dilute domain wall gas, the expression (|4. 1 8fl constitutes an approximate solution to Eq. (|2.6() . Since the 
domain wall solutions are associated with the transient submanifold p = it also follows from Eqs. (|4.3|) . (|3.II|1 and 
(|3.IU|) that they carry vanishing energy, momentum, and action within the canonical phase space approach. 



V. DOMAIN WALL DYNAMICS 



In order to impart dynamical attributes to the domain walls and thus provide solutions to the coupled field equations 
j|4.I(l and 1|4.2(1 we perform a linear stability analysis about the static domain wall solutions. 

A. Dynamics of a single domain wall 

Setting u — Udw + Su and p = pdw + <5p, Pdw = in Eqs. 14.1fl and (|4.2() we obtain to linear order the coupled 
equations 

^ = + 2rfc 2 (l - 3ulJ6u + Sp, (5.1) 

^ = -r0-2T*g(l-3^ w )5p. (5.2) 

Noting that in a Schrodingcr equation analogue the domain wall profile Udw gives rise to a Bargmann type potential 
Eqs. (|5.1() and i|5.2|) arc readily analyzed, see e.g. Ref. j3^|. Expanding Su and dp on the eigenfunctions ^> n of the 
Schrodinger operator 



D = -V 2 + 2fc 2 

according to 



cosh 2 kox 



(5.3) 



5u = J2 U n^ni (5.4) 

n 

5p = ^Pn^n, (5.5) 



the time-dependent expansion coefficients u n and p n are determined by the coupled equations of motion 

^jT = ~T\ n u n +p n , (5.6) 
at 

= TX n p n , (5.7) 

where A„ is the eigenvalue in the eigenvalue equation D^ n — X n ^> n . Equations (|5.6|l and (|5.7() have the same structure 
as the Fourier transformed versions of Eqs. I|4.4|) and l|4.5|l in the case of the linear diffusive modes, and the solutions 
are given by Eqs. (|4.7(l and Ij4.8|l with Tk 2 replaced by TA n . The canonical phase space structure with Tk 2 replaced 
by rA„ is depicted in Fig. |3J 

The spectrum of D in Eq. I|5.3(l is composed of two bound states Wo an d with eigenvalues Ao = and Ai = 3/cq, 
respectively, and a band of phase-shifted plane wave solutions Wfc with eigenvalues X k = k 2 + 4/cq 

( 'ik \ 1 I 2 1 

*0 = -T 7T7 — ' A o = 0, (5.8) 

\ 4 ) cosh kox 

3fcp \ sinh kgx 



*i= Xl = 3kl (5.9) 

V * / cosh kox 

/ 1 \ 1/2 

*fe=f^J e ikx s k {x), A fc =fc 2 +4fc 2 , (5.10) 

k 2 + k?. — 3k% tanh 2 k x + 3ikk Q tanhfco^ 
Skix) = r, r-rj-, „ , ; ■ (5. If) 

y ' (k - ik )(k - 2ik ) y ' 
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The complex space-dependent s- matrix Sk{x) gives rise to a space and phase modulation of the plane wave mode. 
For x — y — oo we have Sk(x) — ► 1, whereas for x — > oo we obtain S).(x) — > exp(i(5 (fc)) exp(i<5i(fc)) where the phase 
shifts So — 2tan _1 (fco/fc) and <5i = 2tan _1 (2fco/&) are associated with the depletion of the band due to the formation 
of the bound states and ^1, respectively. In terms of the expansion coefficients u n and p n the Hamiltonian H, 
momentum II, and action 5* in Eqs. I|4.3[) . (|3.11() . and (|3.10() are given by 

H = lj2Pn(Pn- 2TX nUn), (5.12) 

n 

U = ^2p n / dxw dw V*„ (5.13) 

n J 

S=lj2 I dt PnPn, (5-14) 

expressing the dynamics of the various modes. 



1. Domain wall motion 



The mode VPo with eigenvalue Ao = plays a special role since it is associated with the uniform translation of the 
domain wall. For Ao = the equation of motions i|5.6|) and (|5.7I) take the form duo/dt — po and dpo/dt — with 
solutions uq — p t, pq — est. For the field u associated with the mode \l/o we have u = u^w + "o^o- Inserting an d 
uo = Pot we obtain 

u = u dw + p o-t i—j kg 1 \7u dw cx u dw (x - vt), (5.15) 



describing a domain wall propagating with velocity 



3 



1/2 



v = -P0<r[ W J , * = ±1. (5.16) 

The mode "So is the well-known translation or Goldstone mode restoring the broken translation symmetry associated 
with the localized domain wall mode. In the present canonical formulation the translation mode also implies the 
propagation of the domain wall. 

It is also instructive to consider the noise field p which in the dynamical description corresponds to the noise r\ . In 
the case of a single domain wall we obtain from Eq. (|5.5|l a static configuration. It is, however, clear that the noise 
field must move together with the domain wall configurations and we conclude that terms beyond linear order give 
rise to a renormalization of the noise field profile implying a finite propagation velocity, i.e., 

1 (5 17) 

cosh 2 fc (x — vt) 

The noise field is thus localized at the position of the domain wall corresponding to a noise impulse associated with 
the formation of the domain wall. In Fig. 0] we have depicted a single moving right hand domain wall (index a = 1) 
with associated noise field. The dynamics of a single domain wall follows from Eqs. 15.12fl . (|5.13|) . and l|5.14|l : we 
have 

Eo = \pl (5.18) 



4fc ^ 1/2 



2 1 

Ho = - P o<7 , (5.19) 

So = \Tpl (5.20) 

In terms of the propagation velocity v we have IIo = (4/3)fco« and we can associate an effective mass (4/3)fco with 
the domain wall motion, 

m = h ; (5.21) 



note that the mass vanishes in the limit of a broad domain wall. 
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FIG. 4: In a) we show a right hand domain wall in u moving with velocity v. In b) we depict the associated impulsive noise 
field p 

2. Deformation and extended modes 

The bound state given by Q5.9|l is odd and accounts for the symmetrical deformation of the moving domain 
wall. The time-dependence of u is given by Eqs. 1|4.7|) and l|4.8|l with rfc 2 replaced by TAi = 3rfc 2 . The band of 
plane wave solutions in Eq. 15.10[1 corresponds to corrections to the leading and trailing edges of the domain wall. 
The modes have diffusive character with growing and damped time behavior and are given by Eqs. I|4.7I) and l|4.8|) 
with rfc 2 replaced by T(k 2 + 4k 2 ); note that in the limit fco — > we recover the linear case. The extended modes are 
phase-shifted 2 tan _1 (fco/fc) + 2tan (2fco/k) which by Levinson's theorem corresponds to the two bound states, the 
translation mode and deformation mode, being depleted from the band. 



B. Dynamics of multi-domain wall configurations 



In the multi-domain wall case the analysis proceeds in a similar manner. Expanding about the static multi-domain 
wall configuration in Eqs. (|4.18|) and (|4.19l) . u — u^w + Su and p = pdw + Sp, Pdw — 0, we obtain for a dilute domain 
wall gas the coupled linear equations 



dS u 
dSp 



= r 



Q X 2 + 2LI *0 



-2 + 3£ 



1 



= -r 



^ - 2rfc 2 

dx 2 



-2 + 3E 



^ cosh ko(x — Xi) 
1 



j cosh fco (x — Xi) 



Su + Sp, 
dp, 



which are analyzed in terms of the spectrum of the Schrodinger operator 

1 



D = -V 2 + 2k 2 



2-3E 



j cosh fco (a; — x{) 



(5.22) 
(5.23) 

(5.24) 
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with identical well-separated potential wells at X4. Expanding Su and Sp on the eigenfunctions \&„ of D, D$ n = X n ^> n , 
we recover the equations of motion l|5.6[) and l|5.7|l . 

The eigenstates are readily expressed as linear superpositions of the eigenstates for the individual potential wells 
and we obtain the translation modes 'J'o with eigenvalue Ao = 0, the deformation modes \Pi with eigenvalue Ai = 3k 2 , 
and the extended plane wave modes with eigenvalues A& = k 2 + 4k 2 , 

n 1 

* = A o*l % * 7JT , z, (5-25) 

cosh ko\x — Xi) 



sinh ko(x — Xi) 
3sh 2 k$(x — xi) 

k 2 + k 2 — 3/cq tanh z k$x + 'iikko tanh k^x 



* 1 = $>i*i, ttioc ^rf'^ , (5-26) 

~^ cosh kq{x — xi) 



ihf TT / \ / \ t + Kn - OK n liUlIl KqX -f OIKKQ IHU.IL KQX /v«w> 

^oce^n = (fc-ifaXfc-Mfo) ■ (5 - 2?) 

In terms of the individual eigenfunctions the expansions of Su and Sp take the form 

Su = ^ u ni^n, Uni = U n A % , (5.28) 
III 

Sp = Y2pm^ni u ni =p n Ai, (5.29) 



and we obtain the equations of motion 



dt 

dpni 

dt 



= -TX n U n i+p n i, (5.30) 

= TX„p ni , (5. 31) 



with solutions given by Eqs. (|4.7|l and l|4.8(l . Moreover, the dynamics of a multi-domain wall configuration is given 
by 



H = \Y.P*m(Pm ~ 2rA„u m ), (5.32) 

ni 

U = J2p ni J dxu dw W l n , (5.33) 

ni 

\Y,j dt P*mPm- (5-34) 



S 2 

ni 



For the translation mode in particular we have du^i/dt — poi, dpoi/dt = 0, yielding UQi — po{t and we obtain 

1/2 



U = Itch 



(3k \ 
— J fc^" 1 Vudw oc Udw(x - Vit), (5.35) 



with propagation velocity 



3 ^ 1/2 



-Po°i ( ) • (5-36) 



Likewise, the dynamics of the domain wall at X{ is given by 

1 



Em = (5.37) 

1/2 



5 = \Tp 2 0i . (5.39) 
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FIG. 5: We show the transition from the ground state u = —1 to the ground state u — + 1 in a system of size L in time T due 
to the nucleation of a domain wall pair at the center of the system. 

Summarizing, the linear analysis of the static domain wall configuration leads to a picture of a dilute gas of propagating 
domain walls. Superposed on the domain walls are localized deformation modes and extended modes of diffusive 
character. The time evolution of the domain wall-linear mode gas is moreover subject to three constraints: 1) the 
topological signature of the domain walls implies that a right hand domain wall is matched to a consecutive left hand 
domain wall, 2) translational invariance implies that the total momentum II given by Eq. (|5.33|l is conserved, and 3) 
finally, time translation invariance entails that the total energy E given by Eq. i|5.32[l is a constant of motion. 

C. Domain wall nucleation and annihilation 

The above analysis of the dynamics of the domain wall - diffusive mode system is restricted to the dilute gas regime. 
However, the topological constraints together with the conservation laws allow a heuristic analysis of domain wall 
collisions. When a right hand domain wall collides with a left hand domain wall the topological constraint implies 
that they must annihilate to the uniform ground states u = +1 or u = — 1. Moreover, since Eq. I|3.11fl implies that the 
momentum II = for u = ±1, the domain wall pair prior to collision must move with equal and opposite momenta, 
i.e., equal and opposite velocities. Since the phase space formulation is time reversal invariant we also infer that pairs 
of oppositely moving domain walls of opposite parity can nucleate out of the uniform ground states u = ±1. In the 
ground state Eq. 1)4. 3|) implies that the energy is given by E — (1/2) J dx p 2 . Consequently, in order to generate 
domain wall pairs out of the ground state we must assign a finite noise field p. In Fig. [S]we have in a plot of t versus 
x depicted the transition from the state u = — 1 to the state u — +1 due to the formation of a domain wall pair. The 
system is of size L and the transition takes place in time T. In Fig. [S]we have shown the corresponding propagation 
of the domain wall pair receeding from the nucleation zone with opposite velocities together with the associated noise 
field. According to Eqs. I|5.16[) and l|5.17|l the noise field profiles comoving with the domain walls are positive and 
have the form 1 / cosh 2 fco (x — vt) 

VI. STOCHASTIC INTERPRETATION 

The domain wall gas picture introduced above allow a systematic dynamical approach to the determination of 
kinetic pathways and to the evaluation of the Arrhenius factors associated with the transitions. The approach 
moreover permits a straightforward stochastic interpretation making contact with the customary discussion of the 
noisy Ginzburg-Landau equation. 
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FIG. 6: In a) we show the domain wall pair propagating away from the nucleation zone. In b) we show the associated comoving 
noise field profiles. 
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FIG. 7: We show the orbits in the («o,Po) phase space in the case of domain wall motion. The finite-time orbits lie on the finite 
energy manifolds Eo = (l/2)pg- In the long time limit T — > oo the orbit from ui to 112 migrates to the zero-energy manifold 
E = 0. 

A. Domain wall random walk 

In the case of a single domain wall the dynamics follows from Eqs. I|5.ti[) and (|5.7|l for Ao = 0, i.e., the equations 
of motion dua/dt = po and dp$/dt — with solutions uq = p^t, po = constant. The energy Eq, momentum IIo and 
action So are given by l|5.18|l . (|5.19|) . and 15.20|l . respectively, and the orbits lie on the energy manifolds E = (1/2)j3q. 
In Fig. we have depicted the phase space for the motion of a single domain wall. The plot depicts orbits from m 
to U2 in time T on the manifold Eq = (1/2)pq. In the long time limit the orbits migrate to the zero-energy manifold 
E = 0. The phase space plot is a degenerate limit of the phase space plot in Fig. for ko — > 0. We note the absence 
of a stationary zero-energy manifold and saddle point yielding a stationary state. Since the momentum IIo according 
to Eq. H5.19fl is proportional to the canonical momentum po it follows directly from the canonical structure that the 
conjugate variable to IIq is the center of mass position xq of the domain wall. From the Poisson bracket {uq,Pq} = 1 
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and {zo,n } 
of motion 



with solution 



1 together with IIq 



-cr(4/co/3) 1 / 2 po J we infer u 
dxo 



x = 



dt 



5°, 

m 



n , 



constant. 



am^^Xo, m = 4fco/3, and the ensuing equation 



(6.1) 



Likewise, the energy and action are given by Eq = IlQ/2m and 

TT 2 



(6.2) 



(6.3) 



The stochastic interpretation of the phase space dynamics readily follows. Considering an orbit from x\ to X2 traversed 
in time T, corresponding to the propagation of a domain wall with center of mass position x\ at time t = to the 
center of mass position X2 at time t — T, and inserting the solution l|6.2[l in Eq. (|5.20l) we obtain 



1 (x 2 -x 1 ) 2 
So = ~ 2 m - , 

yielding the transition probability from x\ to x 2 in the time interval T 

m (x2 — x\) 2 



P(xi,X2, T) oc exp 



2A 



T 



(6.4) 



(6.5) 



This is the Gaussian distribution for random walk with root mean square deviation x rms = (2TA/m) 1 ^ 2 , m — 4fco/3. 
We note that the distribution in Eq. I|6.5|) is obtained in the limit T — > from the overdamped oscillator distribution 
in Eq. I|4.13|l . We conclude that the uniform or ballistic motion of a domain wall within the dynamical description 
corresponds precisely to ordinary random walk of the domain wall within the associated stochastic description. 

The random walk behavior of a domain wall also follows easily from the Ginzburg Landau equation 1)2. 1|) . Inserting 
the fluctuating domain wall ansatz u(x, t) = tanh/co^ — x (t))i where x(t) is the time-dependent center of mass x and 
noting that 5F/Su = for a domain wall we obtain, integrating over space, setting rj(t) = J dx rj(x, t), 



dx(t) 
dt 



OC T)(t), 



(6.6) 



which is the Langevin equation for random walk. 

It is also straightforward to include the contribution to the stochastic behavior from the deformation and diffusive 
modes associated with the domain wall. From Eq. (|4.12() applied to the local deformation mode (def) ui and the 
diffusive modes (diff) it& we obtain for the total action for a dressed domain wall 



S — So + Sdcf + Sdis, 
where So is given by Eq. (|6.4() and Sdm and St, by 



Sdcf — 3r/cQ 



(U 2 1 



une 



-3TkiT\2 



1 



s', !iff = / ^T(k 2 



2tt 



4fc 2 ) 



\u2k-u lk e-^ k2 +0 T \ 2 



(6.7) 

(6.8) 
(6.9) 



For the transition probability from an initial dressed domain wall configuration u\{x) — {x\, itn, uifc} to a final 
dressed domain wall configuration 1/2(2;) = {x2, M21, «2fc} during the time interval T we finally obtain 



P{ui,u 2 , T) = Po-Pdcf-Pdiff, 



(6.10) 



where Pq is given by Eq. I|6.5I) and Pdef oc exp[— Sdef/A] and Pdiff oc exp[— Sdig/A]. In the long time limit t — > 00 the 
domain wall performs a random walk, whereas the deformation mode and diffusive modes attain a stationary state, 
i.e., 



P(a;,T;ui,{wfc)}) cx exp 



rax 



2 1 



2AT 



exp 



3Tk 2 ui 



exp 



1 

'A 



dk. 
27 



r(fc 2 + 4fc 2 )| Ufc | 



(6.11) 
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B. Domain wall gas 



In the case of a dilute domain wall gas with associated linear modes the discussion above applies in a generalized 
form. From Eqs. (|5.3U[) and i|5.31|) for Ao = we obtain the equations of motion duoi/dt = poi an d dp oi /dt — with 
solutions u oi — p oi t + est, p oi — est. The energy, momentum, and action are given by Eqs. (|5.37|) . i|5.38[l . and i|5.39[l . 
respectively, and each domain wall lies on the corresponding energy manifold E$i = (1/2)^. For each domain wall 
the phase space plot in Fig.0|thus applies. Since the action in (|5.34|) is additive with a contribution from each domain 
wall the transition probabilities factorize and we obtain for the random walk part from (|6.5(l 



P({x l i}, {x i2 },T) oc ]^[exp 



m(xi2 - x a f 



2AT 



(6.12) 



Likewise, the contributions from the deformation and diffusive modes follow from Eqs. I|tj.8|) arid l|tj.9fl . 



VII. KINETIC TRANSITIONS 



In the previous sections we have established the dynamical framework in the case of the noise-driven Ginzburg- 
Landau equation and established the connection to the customary stochastic interpretation. In summary, the domain 
wall gas picture provides a description of a switching scenario in terms of moving right hand and left hand domain 
walls with associated linear modes. The dynamical approach moreover implies that the kinetic pathway from an 
initial configuration to a final configuration is associated with a specific orbit in canonical phase space. The action 
associated with the orbit yields the Arrhenius factor associated with the transition. The domain wall gas picture 
shows that a class of orbits, i.e., a class of solutions of the field equations can be parametrized in terms of a dilute gas 
of propagating domain walls with superimposed local deformation and extended phase-shifted diffusive modes. This 
picture is derived from a linear analysis and only holds a priori in a dilute gas limit and at short times. As discussed 
below the domain wall dynamics at later times can be extracted heuristically from a combination of selection rules 
and conservation laws. 



A. Equilibrium 

Before we discuss kinetic transitions it is instructive first briefly to discuss the equilibrium properties. In equilibrium 
the free energy landscape is determined by the structure of F(u) in Eq. 1)2.3(1 . The complex landscape is characterized 
by the presence of two global minima corresponding to the uniform ground states u = +1 and u = — 1 and an infinity 
of saddle points corresponding to static nonuniform multi-domain wall configurations. The flat parts of the saddle 
points correspond to the translation modes associated with the domain walls. The free energy F — for the ground 
states. For an n-domain wall configuration F = nm, where the mass m = 4fco/3. In Fig. [H] we have in sketched the 
free energy landscape associated with a single static domain wall connecting u = — 1 and u = +1. In the case of zero 
noise for A = the system starting from an arbitrary initial configuration ui(x) approaches the state with lowest 
free energy compatible with the imposed boundary conditions. The relaxational dynamics is governed by Eq. (|4.15() 
and corresponds in the dynamical approach to an orbit confined to the p = zero-energy submanifold. For open or 
periodic boundary conditions the configurations with lowest free energy are the two degenerate uniform ground states 
u = ±1. Imposing vanishing boundary conditions u — at x — and x = L for a system of size L the field u grows to 
the uniform values u = —1 oru = +1 over a healing length of order l/fco- Imposing periodic boundary conditions we 
infer that the healing profile corresponds to a half domain wall and that the free energy of the configuration equals 
2 x (l/2)m = m. Finally, for skew boundary conditions with u = —1 for x — and u = +1 for x — L the state with 
lowest free energy corresponds to a single domain wall with free energy m. The different scenarios are depicted in 
Fig. [5] In the presence of noise for A ^ the ground states become metastable and the system populates the excited 
states with F ^= 0. The partition function provides a global characterization of the equilibrium excursions in the free 
energy landscape and is according to Eq. 1)2. 5fl given by 



^exp 



2F 



(7.1) 



where the configuration {u} has the free energy F(u), yielding the Boltzmann factor in Eq. I|7.1|l . The evaluation of 
Z and associated correlations, e.g., {uu)(x,t) in the dilute non-overlapping domain w all g as limit follows closely the 
well-known soliton and instanton methods developed in the seventies and used in Ref. |34| for the classical easy-plane 
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FIG. 8: In a) we have in a contour plot sketched the free energy landscape in the case of a single static domain wall overlapping 
for large \x\ with the ground states u = — 1 and it = +1. The points A and B indicate two distinct positions of the domain wall, 
corresponding to the translation mode. In b) we have sketched the free energy as function of u. The u axis labels collectively 
the u(x) configurations in the free energy landscape. The points A and B refer to the domain wall positions. Since the free 
energy of the domain wall is independent of the center of mass position the landscape exhibits a constant F ridge. In c) we 
have depicted the corresponding domain wall configurations at center of mass positions A and B. 



ferromagnet. Expanding the free energy defined by Eqs. I|2.3fl an< i H2.4|) about an n-domain wall configuration we 
have, setting u = iidw + Su, F — F(ud w ) + (l/2)SuD8u. Using F(ud w ) = rim, inserting D from Eq. I|5.24fl . an d 
expanding Su on the eigenfunctions we obtain 



F 



2 ^ 

ni 



K\u r . 



(7.2) 



Introducing a domain wall chemical potential fi^w m order to control the domain wall densities we arrive at the grand 
partition function 



Z = 2 ^2 ex P 

n=0 



2T 



1 



Y[ du 0l du u Y[ duki exp 



r 

A 



^0 U 0i + A l 

\ 1=1 



1=1 



fc 



\Uki 



(7.3) 



The overall factor 2 arises from the double degeneracy of the ground state and the domain walls connecting the ground 
states. The Ao = eigenvalue yields the translation mode and it follows that du^i — m}/ 2 dxi, where Xi is the position 
of the i-th domain wall. The factor 1/n! takes into account the ordering of the domain wall when integrating Xi over 
a system of size L. Performing the Gaussian integrals over the deformation and diffusive modes we have in more 
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+ 1 



FIG. 9: In a) we have shown the ground state configurations with F = in the case of open or periodic boundary conditions. 
In b) the lowest configuration with fixed boundary conditions have F — m. In c) we show the lowest spatially degenerate 
domain wall configuration with F — m in the case of skew boundary conditions. 



reduced form 



n=0 



exp 



2r 



(mA)"/ 2 L" 



vrAA 

Ta7 



i/2 



exp 



1 v-^ , ttAA 

2? 1os Ta7 

k 



(7.4) 



We have introduced the large wavenumber UV cutoff A in order to regularize Z at small distances. Note that in the 
treatment in Ref. [34[ the lattice distance of the magnetic chain provided a natural small scale cutoff. Replacing the 
summation over k by an integral according to the prescription — > J dkpk, where the density of states pk = L/2tt + 
(n/2ir)d6/dk, 6 = 2 tan _1 (/so//j) + 2 tan~ 1 (2fco/fc), and summing over domain walls the partition function factorizes 
into a diffusive part and a domain wall part, incorporating also the contribution from the localized deformation modes, 



ZdiffZdw, 



(7.5) 
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FIG. 10: We show the domain wall density nd w in units of A/2r as function of the noise strength A in units of Tm. The 
density vanishes in the limits A — » and A — * oo and exhibits a maximum at 4mF/3. 



Zdiff = exp 



LA, ttA 

2^ l0g rA 



Z dw = 2 exp 



8Lk, 



3/2 



r\ 1/2 (2T 



(7.6) 
(7.7) 



In a thermal environment A = 2TT and the partition function gives direct access to e.g., the specific heat. We shall 
not pursue this calculation here except noting that the gap in the domain wall excitation spectrum gives rise to a 
Schottky anomaly in the specific heat, see e.g., Ref. (3^|. Noting that the number of domain walls accessed by the 
stationary fluctuations is undetermined it is, however, instructive to evaluate the domain wall density n^w- From the 
structure of Z we infer n dw = (A/2r)(dlog Z/c?^) Al=0 an d inserting Eqs. (|7.5|) . (|7.6[l . and (|7.7|l we obtain, setting 
fcp = im/A,x — > / 



8tt 



-1/2 



3r 

4A' 



3/2 



exp 



2r 



(7. 



The domain wall density vanishes in the limits A — > and A — > oo and exhibits a maximum for A « Tm. We have 
depicted function of A in Fig. ITUl 



B. Transitions in the Kramers case 



Before turning to the noise- induced kinetic transitions in the Ginzburg-Landau equation it is instructive to review the 
classical Kramers theory 0, l37j for a single degree of freedom, which can be analyzed completely, and its formulation 
within the canonical phase space approach. 



Kramers theory 



We consider the overdamped case characterized by the Langevin equation for one degree of freedom u 

du ^dF .„ „. 

-^ = - T -r + v, (7.9) 

at du 

( V (t)v(0)) - AS(t), (7.10) 
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with kinetic coefficient T and noise strength A; in thermal equilibrium A = 2TT, where T is the temperature. For the 
free energy F we assume a double well profile with maximum -F(O) at u — and minima at u = ±1 with F(±l) = 0. 
The free energy is depicted in Fig. El a). The Fokker-Planck equation associated with Eqs. (|7.9|1 and l|7.1t)|l has the 
form 

dP 1 . d 2 P d f„dF 
We note the stationary state for dPj dt = 



Qt=2^ + Tu r ^ P )' (7 ' n) 



exp 



(7.12) 



in accordance with Eq. I|2.5[) . showing how the free energy profile is globally sampled in the stationary state. 

In order to evaluate the transition rate from, say, the state u — — 1 to the state u = +1 across the free energy barrier 
F(0) we follow Kramers and set up a constant probability current J across the barrier generated by a source to the 
left of the barrier in region A and absorbed by a sink to the right of the barrier in region B. From the conservation 
law dP/dt + d J I du = and Eq. (|7.11|l we derive the current 

1 BP d,F 

J = -t A ^--T^-P. (7.13) 

2 du du y ' 

Assuming a sink at u + > we obtain the steady state solution 

P(u) = ^ exp (-^f («)) £ + du' exp (^F(u')j . (7.14) 
Setting up a population n = duP(u) to the left of the barrier the rate is given by k — J/n and we have 



2 r° f 2Y \ r u+ f 2T \ 

fc-^-y duexp(-—F(u)) J du'vq>(—F(u')\. (7.15) 

Finally, a simple steepest descent calculation for A — > yields Kramers celebrated result for the rate in the overdamped 
Smoluchowski limit, F" = d 2 F/du 2 , 

k = ^"(-l^'WD^exp (-f^(O)) • (7.16) 

The rate is determined by the Arrhenius factor exp(— (2r/A)F'(0)) depending on the height of the barrier and the 
prefactor (T/2tt)(F" ( — 1)|F"(0)|) 1 / 2 . Here the double derivative F" can be associated with the oscillation frequencies 
in the potential well at u = — 1 and about the maximum at u = 



2. Dynamical interpretation of Kramers theory 



Here we discuss the Kramers escape problem from a potential well within the canonical phase space method. 
According to the general formulation in Sec. Ill the Hamiltonian associated with the Langevin equation (|7.9l has the 
form 



yielding the equations of motion 



and the associated action 



du _ ^dF 

dt du ' 

dt du 2 ' 



/u,T 
dt 



du 

p Tt~ H 



(7.17) 

(7.18) 
(7.19) 

(7.20) 
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FIG. 11: In a) we show the double well structure of the free energy F. In b) we show orbits in the canonical phase space. The 
zero-energy manifolds are p = and p = 2YdF/du, intersecting at the saddle points (u,p) — (— 1, 0), (0, 0), and (1, 0). In c) we 
show the inverted double well potential entering in the Newton equation description of the transition. 
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For the normalized transition probability we have 



P(u,T) 



exp 




S(u,T)~ 
A 




J ditexp 


S(u,T)~ 
A 



(7.21) 



In Fig. El b) we have depicted the phase space for the Kramers escape case. The zero-energy manifolds are given by 
p = and p = 2TdF / du. On the phase space plot we have shown an orbit from region A to region B, i.e., across the 
free energy barrier. Comparing Fig.^Ja) with Fig.^Jb) we note that the "uphill" part of the orbit in region A is 
controlled by the p = 2TdF/du manifold whereas the "downhill" part in region B is controlled by the p = manifold. 
In the long time limit the orbit approaches the zero-energy manifolds p — and p — 2TdF/du passing close to the 
saddle points at u = — 1, u = 0, and u = 1. A simple calculation along the "bulge" setting H — and inserting 
p = 2TdF/du in Eq. (|7.2U|) yields the Arrhenius factor 



exp 



2r 



F(0) 



(7.22) 



in accordance with Kramers result in Eq. H7.16(l . 

Finally, eliminating p in the coupled equations l|7.18|l and (|7.19|) we obtain the Newton equation of motion 



a u 

lit 2 



du 



1 fdF 
2\~du 



(7.23) 



for the motion of a particle of mass 1/T 2 in the potential — (1 / 2) (dF / du) 2 . In Fig.^Jc) we have shown the potential 
which for a double well free energy possesses three maxima at u — and u = ±1. As indicated in Fig. II II c^) the long 
time orbit from A to B is associated with the long waiting time at the maximum for u = 



C. Transitions in the Ginzburg-Landau case 

In order to induce a non-equilibrium kinetic transition across the free energy landscape in the Ginzburg-Landau 
case we fix the initial configuration ui(x) at time t = and the final configuration 7/2(2;) at time T. There are 
two fundamental issues: 1) the determination of the kinetic pathways and 2) the evaluation of the transition rate. 
Generally, in order to minimize the free-energy cost, the pathway passes via saddle points, i.e., the multi-domain wall 
excitations, in the free energy landscape. The actual path chosen, however, depends on the time T allocated to the 
transition. The transition rate is determined by the Arrhenius factor exp(— 2rF/A) associated with the free energy 
F of the saddle points encountered along the pathway. In addition there is a prefactor determined by the attempt 
frequencies; this issue, however, is not dealt with in the present context. 

Focussing on a kinetic transition or switch between the two ground states u = — 1 and u = +1 the pathway 
corresponds to nucleation of domain walls and their subsequent propagation across the system. Let us first consider 
the case of fixed boundary conditions u — at x — and x = L. The ground state configurations are depicted in 
Fig. 0b). In the presence of noise the ground states are metastable and the system can switch from u = — 1 to u = +1 
by means of the propagation of a single domain wall across the sample in the transition time T. As discussed above 
the profiles close to the boundaries over a healing length of order fc^ 1 correspond to half domain walls and the free 
energy of either ground state is equal to Fa w = m, m = 4fco/3; note that for periodic boundary conditions as depicted 
in Fig. a) the free energy of the ground states vanishes. In order to effectuate a transition a domain wall must 
nucleate at either boundary and subsequently propagate across the system. The nucleation free energy is F^w = fn 
and from our discussion of the Kramers case we infer the Arrhenius factor exp(— 2rm/A) or within the dynamical 
description a nucleation action of the order 

5 nucl = 2rm. (7.24) 

Within the conventional stochastic description the noise generates both the nucleation and the subsequent diffusion 
of the domain wall across the sample. Since the domain wall carries a finite energy we must within the dynamical 
description assign a finite energy or finite noise field p(x) at t = in order to ensure propagation of the domain wall. 
Considering the case where the right hand domain wall nucleates at x = 0, we assign the noise field p{x) oc po/ cosh 2 k$x 
given by Eq. H5.17[) at time t = 0. In terms of the momentum LTo = — pom 1 / 2 given by Eq. (|5.19|) we thus obtain the 
velocity v = U /m, energy E = IIo/2m, and action by Eq. I|tj.3|) . i.e., Sq = TUq/2iti. The constraint that the system 
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FIG. 12: In a) the switching from u = +ltou=— lin time T is effectuated by means of a right hand domain wall propagating 
with velocity v = 1/T. The domain wall is nucleated at x — and annihilated at x = 1. In b) the process is depicted in an 
(x,t) plot. 



of size L is switched at time T moreover impose a selection rule on the propagation velocity. We infer v = L/T and 
obtain the action So = mL 2 /2T associated with the transition. The total action associated with a single domain wall 
switch is thus given by 



5*1 (T) — Snucl + 



mi 
2T 



(7.25) 



In Fig. ll2l we have in a) depicted the propagation of a single right hand domain wall across the sample. In b) we have 
shown the transition in an (%, t) plot. In the present case with fixed boundary conditions the momentum IIo associated 
with the motion of the domain wall is generated at the boundary x = at time t — and, subsequently absorbed at 
the boundary x — L at time t = T. The momentum is given by Eq. (|5.13|) . i.e., n = po J dxudwd^o/ dx. Inserting 
and «dw given by Eqs. I|5.8[l and (|2.7|l . respectively, and performing a partial integration using J dx cosh -4 fco^ = 4/3fco 
we obtain IIo = —pom 1 / 2 > since po < for a positive propagation velocity Correspondingly, the momentum is 
absorbed at x = L. 

For fixed boundary conditions the transition from u = — ltow = +l can also take place by means of nucleation 
of two domain walls at both ends of the system. The domain walls subsequently propagate towards one another and 
annihilate at x — L/2. The scenario is shown Fig. 1131 Since the action is additive for the two-domain wall system the 
nucleation action is given by 2S' lluc i = 2m. However, in order to effectuate the transition in time T the velocity is half 
the velocity in the single domain wall case and we obtain for the action associated with the 2-domain wall transition 



£2 {T) — 2Snucl " 



ml? 
AT 



(7.26) 
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FIG. 13: In a) the switching from u — +1 to u = —1 in time T takes place by means of two domain walls propagating in 
opposite directions with velocity v — 1/2T. The domain walls nucleate at the boundaries and annihilate at the center. In b) 
the switching process is depicted in an (x,t) plot. 



Finally, in the general case of a transition from m = +1 to ti = - lin time T by means of the nucleation and 
propagation of n domain walls the nucleation action is n5 nuc i and the domain walls move with velocity L/Tn; we 
obtain the action 

S n {T) = nS nncl + ^. (7.27) 

In Fig. 1141 we have in an (xt) plot shown a 4-domain wall switch. 

In addition to domain wall modes time-dependent localized deformation and extended diffusive modes are also 
excited, corresponding to small Gaussian fluctuations about the local minima and saddle points and the transition 
pathway from u = +1 ton = —1 proceeds by propagating domain walls with superposed linear modes subject to energy 
and momentum conservation and topological constraints. The energy of the initial state is given by E = (1/2) J dxp 2 
and the noise field thus has to be assigned initially in order to reach the switched state u — —1 in a prescribed time T. 
For topological reasons the domain walls must nucleate and annihilate in pairs subject to absorption or radiation of 
linear modes, respectively. Since the linear modes also carry positive action the dynamical modes with lowest action 
correspond to nucleation or annihilation of domain wall pairs with equal and opposite momenta, i.e., equal velocities. 

In the case of periodic boundary conditions the momentum II = J dxudp/dx of the initial and final states is zero. 
The system is translational invariant and the formation and annihilation of one or several domain wall pairs moving 
with the same speed take place at equidistant positions along the axis. For fixed boundary conditions the translational 
invariance is broken and the momentum II is nonvanishing corresponding to nucleation and annihilation of domain 
walls at the boundaries. This general scenario of switching is completely consistent with the numerical analysis in 
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FIG. 14: We show a 4-domain wall switch in an (x,t) plot. 




1. Improved estimate of S nuc i 



The above estimate of the nucleation action S nU ci for a domain wall was based on an analogy with the Arrhenius 
factor in the simple Kramers case of escape from a potential well. Here we improve the estimate of S nuc i based on 
the phase space formulation. We consider the case of a switch from u = +1 to u = —1 in time T proceeding by 
1) the nucleation of two domain walls at the center of the system during the time interval 5t, 2) propagation of the 
domain walls with opposite velocities during the time interval T — 26t, and 3) the annihilation of the domain walls at 
the boundaries during the time interval St. From the width of a single domain wall k^ 1 and the propagation velocity 
v = L/2T we estimate the nucleation time, i.e., the time it takes for the nucleation to separate in two distinct domain 
walls, to be of order St w 2/vk = AT/k L. From Eq. I|3.10fl the nucleation action per domain wall is given by 



1 



u.St 



S'nuci = - / dxdt(pdu/dt - H), (7.28) 

where u is the nucleation configuration just prior to breaking up into two domain walls; TL is the energy density 
which we determine below. In order to estimate S'nuci we consider the field equations l|4.1|l and (|4.2|) . The initial 
configuration is u(x) = u(x,t = 0) = 1 and denoting the initial noise field by po(x) = p(x,t = 0) and, moreover, 
considering a large system so that St is small, we obtain to leading order in St from the equations of motion 

,d 2 u(x,t) 
dx 2 



u{x,t) w r v ' ' +2Tk^u(x,t)(l-u(x,tY) +po(x))St + l, (7.29) 



t=o 



p(x, t) « {-T^^- - 2Tk*p(x, t)(l - 3u(x, t) 2 )^j St + p (x), (7.30) 

or by insertion 

u(x,t) &p (x)St + 1, (7.31) 

p(x, ~ r (-0 + 4fc oV>(£)) St + p {x). (7.32) 

These solutions describe the initial part of the orbit in (up) phase space corresponding to the nucleation of two domain 
walls. The unspecified initial noise field Po(x) acts as an initiator. The noise profile Po(x) is localized at the position 
of the nucleation. Since the orbit lies on an energy surface we have Eq = (1/2) J dxpo(x) 2 and we infer the energy 
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FIG. 15: We sketch the orbit in phase space corresponding to the nucleation of a domain wall pair during time t nuc i, the 
subsequent propagation across the system during time T — 2t nuc i , and the final annihilation of the individual domain wall at 
the boundaries during time tnuci- 



density TL = (l/2)po(x) 2 . From the equation of motion du(x)/dt « Po(x) and inserting in the action in Eq. I|7.28|l we 
obtain S^ucl ~ (l/2)<5i / dx\po(x) + T(—d 2 po(x)/dx 2 + Ak' 2 po{x))5t)pQ{x) — (l/2)po(x) 2 ]. Rearranging and performing 
a partial integration, assuming po(x) localized, we find the following expression for the nucleation action per domain 
wall: 

5 nucl « \stEo + ±-(5t) 2 (sklTEo +T [ dx ) ] . (7.33) 



This is a short time estimate and holds for a large system, i.e., for k$L ^> 1. The energy Eq = 2(l/2)mu 2 is found 
from the 2-domain wall sector. Inserting St = 4T/koL, v = L/2T, m — 4fco/3, estimating (dpo(x) / dx) 2 — fk 2 po(x) 2 , 
where / ia a "fudge" factor of order one, and lumping the term linear in i5 with the domain wall pair propagation, we 
obtain 5 nU ci = <5t 2 (4 + f)k 2 TEo. Further reduction yields the result 

S„uci = (4 + f)Tm (7.34) 

This expression for the nucleation action has the same form as the one derived from the simple Kramers theory. In 
Fig. I15l we have in a phase space plot sketched the orbit and indicated the nucleation, propagation, and annihilation 
regions. 



VIII. INTERPRETATION OF NUMERICAL RESULT 



Here we make contact with the numerical analysis of the Ginzburg-Landau equation by E, Ren, and Vanden- 
Eijnden |2l| . These authors analyze the noise induced switching by means of an optimization techniques applied to 
the Freidlin-Wentzel action 

s ™ - s r dx [ dt (s - r £ + 2r *»" (i - ° 2) ) 2 - (tu » 

for a system of size L = 1 over a time span T and find that the global minimum of ^fw corresponds to the nucleation 
and propagation of domain walls. First we observe that the minimizing configurations, the minimizers, of (|8.1I) are 
identical to the orbits found within the canonical phase space approach. This correspondence was discussed in Sec. 
Ill and is seen by noting that insertion of the equation of motion for p and the Hamiltonian (|4.3|l in the expression 
for the action l|3.10|l yields the Freidlin-Wentzel form in (|8.1fl . 
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Switching a system of size L in time T by means of a single domain wall, corresponding to the pathway via the 
lowest local minimum of the free energy at F^ w = m, the propagation velocity v = Po/m = L/T and we obtain the 
action Si(T) — S nuc i + mL 2 /2T in Eq. Ij7.25|l and associated transition probability 

p — (-¥)=>• (-55) ■ <8 ' 2) 

In the thermodynamic limit L — > oo, P — > as a result of the broken symmetry in the double well potential. At 
long times the action falls off as 1/T. At intermediate times t and positions x we have P oc exp(— mx 2 /2At) and 
the domain wall in the stochastic interpretation perform a random walk with mean square displacement 2At/m, 
corresponding to diffusive behavior. In Fig. El a we have shown a domain wall nucleating at the left boundary and 
propagating with constant velocity v = 1/T to the right boundary, where it annihilates. We have used the same 
parameter values as in |2l|. i.e., 5 = T = .03, 2FkQ = <J -1 , T = 7, and a system size L = 1. In Fig. ^]b we have 
plotted the trajectory of the domain wall in space and time. 

The switching can also take place by nucleating two domain walls at the boundaries. These then move at half the 
velocity v/2 and subsequently annihilate at the center. This process corresponds to the pathway via the local saddle 
point of the free energy at Fd w = 2m, and the action is given by 52 (T) = 25 nuc i0 + 25i(4T) in Eq. (|7.26l) . The 
associated transition probability is 

'—(-^M-is)- < 83 > 

Snapshots of this process are shown in Fig. 1131 a and the corresponding space-time plot in Fig. 1131 b. 

Combining the contributions from nucleation and the subsequent domain wall propagation we obtain the heuristic 
expression in Eq. I|7.27[) . where 5 nuc is the action for nucleating a single domain wall and n is the number of walls. 
From the improved estimate of 5 nuc i in Eq. (|7.34|) with fudge factor / w 1 we have 

5 nuc ~ 5rfc . (8.4) 

In Fig. E| we have plotted 5 versus T for n — 1 — 6 domain walls using the parameter values in j2^]. Choosing 5 nuc 
according to jHUjl we find excellent agreement with the numerical results. 

As also discussed in [2l| we note that the switching scenario depends on T. At shorter switching times it becomes 
more favorable to nucleate more domain walls. In the present formulation this feature is associated with the finite 
nucleation or annihilation action 5 nuc . This is evidently a finite size effect in the sense that the action at a fixed T 
diverges in the thermodynamic limit L — > oo, corresponding to the broken symmetry. 

IX. SUMMARY AND CONCLUSION 

In the present paper we have implemented the dynamical phase space approach developed previously for the noisy 
Burgers equation to the case of the one dimensional noise-driven Ginzburg-Landau equation. Based on a linear 
analysis of the static domain wall solutions in the noiseless case we find that the kinetic transitions take place by 
means of propagating multi-domain wall configurations with superimposed local deformation and extended diffusive 
modes. The approach also allows for a determination of the Arrhenius factor associated with the transitions. The 
motivation for the present work is a recent numerical optimization study by E et. al. pH of the noisy one dimensional 
Ginzburg-Landau equation based on the Freidlin-Wentzel theory of large deviations. We find excellent agreement 
both qualitatively and quantitatively with the numerical finding of E. et al. plj . The dynamical approach offers in 
the nonperturbative weak noise or low temperature limit an alternative way of determining dynamical pathways and 
the Arrhenius part of the associated transition rates. 
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FIG. 16: The action S(T) given by (|7.27pi is plotted as a lunction of T for transition pathways involving up to n — 6 domain 
walls. The lowest action and thus the most probable transition is associated with an increasing number of domain walls at 
shorter times, indicated by the heavy limiting curve. The curves correspond to choosing Snuc = 5Tko- 
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